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Characterization of quantum dynamics is a fundamental problem in quantum physics and quantum infor- 
mation science. Several methods are known which achieve this goal, namely Standard Quantum Process To- 
mography (SQPT), Ancilla- Assisted Process Tomography (AAPT), and the recently proposed scheme of Direct 
Characterization of Quantum Dynamics (DCQD). Here, we review these schemes and analyze them with respect 
to some of the physical resources they require. Although a reliable figure-of-merit for process characterization 
is not yet available, our analysis can provide a benchmark which is necessary for choosing the scheme that 
is the most appropriate in a given situation, with given resources. As a result, we conclude that for quantum 
systems where two-body interactions are not naturally available, SQPT is the most efficient scheme. However, 
for quantum systems with controllable two-body interactions, the DCQD scheme is more efficient than other 
known QPT schemes in terms of the total number of required elementary quantum operations. 

PACS numbers: 03.65.Wj 



I. INTRODUCTION 

Characterization of quantum dynamical systems is a cen- 
tral task in quantum control and quantum information pro- 
cessing. Knowledge of the state of a quantum system is in- 
dispensable in identification/verification of experimental out- 
comes. Quantum state tomography has been developed as a 
general scheme to accomplish this task jit]. In this method 
an arbitrary and unknown quantum state can be estimated by 
measuring the expectation values of a set of observables on an 
ensemble of identical quantum systems prepared in the same 
initial state. Identification of an unknown quantum process 
acting on a quantum system is another vital task in coherent 
control of the dynamics. This task is especially crucial in ver- 
ifying the performance of a quantum device in the presence 
of decoherence. In general, procedures for characterization of 
quantum dynamical maps are known as quantum process to- 
mography (OPT) — for a review of quantum tomography see 
Refs. SSIl. 

There are two types of methods for characterization of 
quantum dynamics: direct and indirect. In indirect meth- 
ods, information about the underlying quantum process is 
mapped onto the state of some probe quantum system(s), 
and the process is reconstructed via quantum state tomog- 
raphy on the output states. We call these methods indirect 
since they require quantum state tomography in order to re- 
construct a quantum process. A further unavoidable step in 
indirect methods is the application of an inversion map on 
the final output data. Standard Quantum Process Tomogra- 
phy (SQPT) m Jsl B| and Ancilla- Assisted Process Tomog- 
raphy (AAPT) 13, 11 H belong to this class. On the 
other hand, in direct methods each experimental outcome di- 
rectly provides information about properties of the underly- 
ing dynamics, without the need for state tomography. In the 
last decade, there has been a growing interest in the devel- 
opment of such direct methods for obtaining specific infor- 



mation about the states and dynamics of quantum systems, 
such as estimation of general functions of a quantum state 
ifTHl . detection of quantum entanglement lfl2ll . measurement 
of nonlinear properties of bipartite quantum states jTs}], esti- 
mation of the average fidelity of a quantum gate or process 
1 14, [Tsll . and universal source coding and data compression 
I idl . The method of Direct Characterization of Quantum Dy- 
namics (DCQD) 113, S [H, IMl is the first scheme which 
provides a full characterization of (closed or open) quantum 
systems without performing any state tomography. In this 
method each probe system and the coiTesponding measure- 
ments are devised in such a way that the final probability dis- 
tributions of the outcomes become more directly related to 
specific classes of the elements of the dynamics. A complete 
set of probe states can then be utilized to fully characterize 
the unknown quantum dynamical map. The preparation of 
the probe systems and the measurement schemes are based 
on quantum error-detection techniques. By construction, this 
error-detection based measurement allows for direct estima- 
tion of quantum dynamics such that the need for a complete 
inversion of final results does not arise. Moreover, by con- 
struction, DCQD can be efficiently applied to partial char- 
acterization of quantum dynamics. For example, as demon- 
strated in Refs. in mi], the DCQD scheme can be used for 
Hamiltonian identification, and also for simultaneous determi- 
nation of the relaxation time Ti and the dephasing time T2 in 
two-level systems. A proof-of-principle optical realization of 
DCQD via a Hong-Ou-Mandel interferometer has also been 
reported ll20ll . Recently, direct approaches for efficient par- 
tial/selective estimation of quantum proc esses based on ran- 
dom sampling have been introduced IZM . Application of the 
direct QPT methods to the efficient parameter estimation of 
many-body quantum Hamiltonian systems is also of special 
interest for practical purposes, and will be addressed in an- 
other publication |23i]. 

In this work, we review all known methods for complete 
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characterization of quantum dynamics, and analyze the re- 
quired physical resources that arise in preparation and quan- 
tum measurements. To the best of our knowledge, this is 
the first complexity analysis of different QPT schemes. We 
conclude that, for quantum systems with controllable single- 
and two-body interactions, the DCQD scheme is more effi- 
cient than the other known QPT schemes, in the sense that 
it requires a smaller total number of experimental configura- 
tions and/or elementary quantum operations. However, for 
quantum systems where two-body interactions are not natu- 
rally available (e.g., photons), the DCQD scheme and (non- 
separable) AAPT cannot be implemented or simulated with 
high efficiency, and the SQPT scheme is in this case the most 
efficient. 

The structure of this paper is as follows. In Sec. HIl we 
briefly review the concept of a quantum dynamical map. In 
the subsequent sections, |llll |IV] and |V] we provide a review 
of the SQPT, AAPT, and DCQD schemes, respectively. Since 
SQPT has been extensively described in earlier literature, we 
provide more detail about the AAPT and DCQD schemes. 
Specifically, we provide a comprehensive discussion of the 
different alternative AAPT measurement strategies, i.e., those 
utilizing either joint separable measurements, mutually unbi- 
ased bases measurements, or generalized measurements. For 
simplicity, we assume that all quantum operations, including 
preparations and measurements, are ideal; i.e., we do not con- 
sider the effect of decoherence during the implementation of 
a QPT scheme. In the final section of the paper — Sec.[yl] — 
we present a detailed discussion and comparison of the differ- 
ent QPT strategies. 



II. QUANTUM DYNAMICAL MAPS 

Under rather general conditions (but assuming a factorized 
initial system-bath state) the dynamics of an open quantum 
system can be described by a completely-positive linear map, 
as follows; 



FIG. 1: Schematic of SQPT. An ensemble of states {pj} are pre- 
pared and each of them is subjected to the map E, and then to the 
measurements \Ern\- 



to the {£";} basis, i.e., characterization of E is equivalent to 
a determination of the independent matrix elements of x, 
where the Ei play the role of observables. When the map E is 
trace-preserving, i.e., — the corresponding super- 

operator X has only — (P independent elements. Hereafter, 
we restrict our attention only to the n-qubit case, i.e., d = 2". 



III. STANDARD QUANTUM PROCESS TOMOGRAPHY 

The central idea of SQPT is to prepare d^ linearly- 
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independent inputs {pk}k^o ^i^^ then measure the output 
states £{pk) by using quantum state tomography Uj, i5|, i60. 
SQPT has been experimentally demonstrated in liquid-state 
NMR fl^j73J^, optical ||23,|2l], atomic IH] and solid-state 
systems [30J. Since the map £ is linear, it can in principle be 
reconstructed from the measured data by a proper inversion. 
Let {pk}k=a ^ linearly independent basis set of operators 
for the space of dx d linear operators. A convenient choice is 
Pk = |m)(n|, where {|m)}j5^Q is an orthonormal basis for 
The coherence |to) (n| can be reconstructed from four popula- 
tions: |m)(n| = |+)(+| + h)(-|-[|m)(7n| + |n)(?7,|](l+i)/2, 
where |+) = (|m) + |n))/v^and |-) = (|m) +i\n))/^/2. 
Linearity of S then implies that measurement of 
f (|— )(— I), £(|m)(m|), and £{\n){n\) suffices for the deter- 
mination of £{\m){n\). In addition, every £{pk) can be ex- 
pressed in terms of a linear combination of basis states, as 
£{pk) = J2i ^kiPi- The parameters X^i contain the measure- 
ment results, and can be understood as the expectation values 
of the fixed-basis operators Ek '■ 



(1) 



Xki = Tr{Ek£ipi)), 



(4) 



where p is the initial state of the system [p E B{H), the space 
of Unear operators acting on TC] and ^l^i =^ I guaran- 
tees that Tt£{p) ^ 1 El. Suppose that {-Bjflo ^ is a set of 
fixed Hermitian basis operators for B{H), which satisfy the 
orthogonality condition 



TriEjE, 



dSi^ 



(2) 



For example, for a multi-qubit system the EiS, can be tensor 
products of identity and Pauli matrices. The Ai operators can 
be decomposed as Ai = "Y^^ a.i,nE„i, and therefore we have 



mn—0 



(3) 



where Xmn = ^miO*nj- The positive superoperator x en- 
compasses all the information about the map £ with respect 



when Ek = Pk- This choice of the Ek is natural, since the pk 
are Hermitian operators and thus they are valid observables. If 
we combine this with the relation EmPkE^i ~ J2i Bmn.ikPi, 
the following equation can be obtained: J2„in Bmn.ikXmn = 
Xki ■ This in turn can be written in the following matrix form: 



Bx = \ 



(5) 



where the x c?^ -dimensional matrix B is determined by the 
choice of bases {pk} and {E„i}, and the d* -dimensional vec- 
tor A is determined from the state tomography experiments. 
The superoperator x can thus be determined by inversion of 
Eq. (|5]), but in general x is not uniquely determined by this 
equation. 

Figure [T] illustrates the SQPT scheme. Let us determine 
the resources this scheme requires. In general, SQPT involves 
preparation of d^ linearly independent inputs {pi}, each of 
which is subjected to the quantum process £, followed by 
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quantum state tomography on the corresponding outputs. As 
we saw above, for each pi we must measure the expectation 
values of the d? fixed-basis operators {Ek] in the state £{pi)- 
Thus the total number of required measurements amounts to 
d^. Since measurement of an expectation value cannot be 
done on a single copy of a system, throughout this paper, 
whenever we use the term "measurement" we implicitly mean 
measurement on an ensemble of identically prepared quantum 
systems corresponding to a given experimental setting. 



IV. ANCILLA-ASSISTED PROCESS TOMOGRAPHY 

In principle, there is an intrinsic analogy between quantum 
state tomography schemes and QPT. This analogy is based 
upon the well-known Choi-Jamiolkowski isomorphism Isill . 
which establishes a correspondence between completely- 
positive quantum maps (or operations) and quantum states, 
£ ^ P£, as follows; 

= (£®/)(|$+)(<I>+|), (6) 

where |$+) = Yl'i=i ^1*) ® I*) maximally entangled 

state of the system and an ancilla with the same size. This one- 
to-one map enables all of the theorems about quantum opera- 
tions directly to be derived from those of quantum states ll32ll . 
In this way, one can consider a quantum process as a quantum 
state (in a larger Hilbert space). ' Therefore, the identifica- 
tion of the original map £ is equivalent to the characterization 
of the corresponding state p£. In other words, the problem 
of quantum process tomography can naturally be reduced to 
the problem of quantum state tomography, and hence, all state 
identification techniques can be applied to the characterization 
of quantum processes as well. The AAPT scheme was built 
exactly upon this basis. 

Generally, within the AAPT scheme, we attach an auxiliary 
system (ancilla), B, to our principal system. A, and prepare 
the combined system in a single state such that complete infor- 
mation about the dynamics can be imprinted on the final state 
{s, 9]. Then by performing quantum state tomography in the 
extended Hilbert space of Hab, one can extract complete in- 
formation about the unknown map acting on the principal sys- 
tem. In principle, the input state of the system and ancilla can 
be prepared in either an entangled mixed state (entanglement- 
assisted) or a separable mixed state. Intuitively, the input state 
in AAPT must be faithful enough to the map £ such that by 
quantum state tomography on the outputs one can identify £ 
completely and unambiguously ifioll . This faithfulness con- 
dition can formalized. Indeed, it is easy to show that a state 
p can be used as input for AAPT iff p has maximal Schmidt 
number, i.e. Sch(p) = i^. ^ The faithfulness condition 



' Choosing {Em} = {\i){j\} results in: x = dpe Si- 

^ Any operator Q acting on a bipartite system AB can be decomposed as 

Q = Ei^ll'^' siAi Bi, where the si are all non-negative numbers, and 
{Ai } and {_B; } are orthonormal operator bases for the systems A and B, 
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FIG. 2: Schematic diagram of separable AAPT. An ensemble of sys- 
tems is prepared in the same quantum state p. Next, they are sub- 
jected to the map £ ® I. Finally the operators {Ej} are measured 
on both the system and the ancilla, which results in the required joint 
probability distributions or expectation values. 



is nothing but an invertibility condition. That is, because of 
linearity of the map £ (g) I, the information is imprinted on 
the elements of the final output states linearly. By performing 
state tomography on the output state {£ (g) I){p), we obtain a 
set of linear relations among possible measurement outcomes 
and the elements of £ and p. The matrix p must be chosen 
such that an inversion becomes possible; thus one can solve 
the set of linear equations for £ I8fl. We provide more details 
below. 

It should be noted that the faithfulness condition is different 
from entanglement. In fact, almost all states of the combined 
system AB (excluding product states) may be used for AAPT, 
because the set of states with Schmidt number less than (P is 
of zero measure. This means that entanglement is not a nec- 
essary property of the input state p in AAPT. Indeed, many of 
the viable input states are not entangled, such as Werner states 
However, it has been argued and also experimentally veri- 
fied that use of maximally entangled pure states offers the best 
performance. That is, even though in principle any faithful 
state can be used in AAPT, the propagation of experimental 
errors from the measurement outcomes to actual estimation of 
{£ /) (p), due to the inversion process, dictates that different 
faithful input states can produce very different errors. 

In order to develop a good faithfulness measure, one can 
consider a general property of a typical input state p, such that 
the output states generated by different quantum dynamical 
processes have maximum distance; e.g., different output states 
should be nearly orthogonal. More specifically, we note that 
the experimental error amplification is related to the inversion, 
which in turn depends on the multiplication by the inverse of 
the eigenvalues of p, s^^- Then, the smaller the eigenvalues 
the higher the amplification of experimental errors. This fact 
has led to the following definition: 

Fip)=Trip^)=J2'i^ 
1=1 

as a proper measure of faithfulness (Ht]. This, indeed, is ex- 
actly the purity of the state p. As a consequence, this implies 
that the optimal (in the sense of minimal experimental errors. 



respectively (34ll . Sch(Q) is defined as the number of terms in the Schmidt 
decomposition of Q. 
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as explained above) faithful input states are pure states with 
maximal Schmidt number and s/ = l/Vd, i.e., maximally 
entangled pure states. 

The Hilbert space of the input state in AAPT is Hab = 
TCa'8)'Hb- At the output, one can realize the required quan- 
tum state tomography by either separable measurements (sep- 
arable AAPT), i.e., joint measurement of tensor product op- 
erators, or collective measurements on both the system and 
ancilla (non-separable AAPT). Both of these measurements 
are performed on the same Hilbert space Hab- Furthermore, 
it is possible to perform a generalized measurement or POVM 
by going to a larger Hilbert space. In the subsequent sections, 
we discuss all of these alternative strategies and argue that 
the non-separable measurement schemes (whether in the same 
Hilbert space or in a larger one) have hardly any practical rele- 
vance in the context of QPT, because they require many-body 
interactions which are experimentally unavailable. 



A. Joint separable measurements 

Let us assume that the initial state of the system and an- 
cilla is pAB - P^JEf ® Ef, where {E^} ({E^}) is the 
operator basis for the linear operators acting on Ha CHb), as 
defined earlier The output state, after applying the unknown 
map £ on the principal system, is the following: 

p'ab = {Sa<» Ib){pab) 

= Y.'^koEk®Ef. (7) 

kj 

In the last line, we have used a^j = J2mniXmnPija'^''''''^ , 
where a™'*'" is defined via E^E^E^^ = J^k 
and depends only on the choice of operator basis. From the 
above equation it is clear that if we consider the basis op- 
erators as observables,'' then the parameters ctkj, which are 
related to the Xnm's, can be obtained by joint measurement 
of the observables E^ (X) E^ . In fact, the expectation values 

Tr{p'j^gE^^ (g) E^^), as the measurement results, are exactly 
the ctkj parameters: 

ak,=Tr{p'^BEf^Ep). (8) 
Now, by defining 

Xki = 2^a^." Xmn, (9) 
rnn 

and considering that the a parameters are known from the 
choice of operator basis, we see that by knowledge of the 



Q!fcj s we can obtain the Xmn parameters through the follow- 
ing matrix equation: 

a = XQ, (10) 

where g [ptj], x = [Xmn], and a = [5^;]. This equation 
implies that unambiguous and unique determination of the x 
matrix is possible iff the g matrix is invertible. After obtaining 
X, by using the linear relation of Eq. (|9]l between x and x 
matrices, one can easily find x by an inversion. 

Equation ( fTOl i implies that if we were to choose p as a mul- 
tiple of the (P X d? identity matrix /, then the unknown quan- 
tum operation, x, would be directly related to the measure- 
ment results, a, without the need for inversion. However, 
positivity of the density matrix pab disallows this choice. For 
example, in the qubit case, it can be easily seen that the opera- 
tors g^ \ l results in pab = \{Ia®Ib+ Xa i^Xb + Ya<5Si 
Yb + Za (8 Zb), which is physically unacceptable because of 
its negativity. Conversely, Eq. ( fTOl i implies that in AAPT no 
choice of the initial density matrix pab can result in a direct 
(inversion-free) relation between the measurement results and 
elements of the unknown map. 

Next, we explicitly show that the invertibility condition of g 
in Eq. ( fTOl i is equivalent to the condition of maximal Schmidt 
number in the corresponding pab- In general, an operator 
Qab on Ha ^Hb can be written as Q = J2jk QjkCj ® Dk, 
where {Cj} {{Dk}) is a fixed orthonormal basis for the space 
of linear operators acting on Ha (Hb)- A singular value de- 
composition of the matrix Q = [Qjk] yields Q = USV, 
where U and V are unitary matrices and is a diagonal ma- 
trix with non-negative entries Sjk ~ SkSjk {sk are the singular 
values of the matrix Q). Using this decomposition, we have 
Q ^ ^i^i Bi, where the operators {Ai = Y,j UjiCj} 
and {Bi = J^k ^ikDk} are also orthonormal bases. This is 
the Schmidt decomposition of the operator Q ifsill . In our 
case, Q is the matrix g. We know that g is invertible iff 
none of its singular values is zero, i.e. VZ, s/ ^ 0. This, in 
turn, guarantees that in the Schmidt decomposition of pab 
(counterpart of Q) all terms are present, that is, it has max- 
imal Schmidt number This confirms that the invertibility 
condition — which is necessary for the applicability of input 
states in AAPT — is exactly what was already termed faith- 
fulness above (for more detail see Ref. |19;]). In fact, even 
separable Werner states, ~ ^^-^ + )('^' I (in which 
\^-)ab = (|01)-|10))Ai3/V2)fores; ^iH, have max- 
imal Schmidt number. Therefore even classical correlation 
between the system and the ancilla is sufficient for AAPT. 

B. Mutually unbiased bases measurements 

Ancilla-assisted quantum process tomography can also be 
performed by using "mutually unbiased bases" (MUB) mea- 



^ If dim('H) = d, then one can choose Eq = I (d X d identity matrix) 
and Ej (j = I, . . . , d^ — I) to be traceless Hermitian matrices. 



Eq. (9) can be written formally as: x = X-^- where x (x) is a row matrix 
obtained by an'anging elements of x (x) in some agreed order, and A is a 
matrix obtained by the corresponding reordering of the fj™''-" parameters. 
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TABLE I: A partitioning of the 2-qubit Pauli group such that the 
eigenvectors constitute a MUB. 



MUB 1 






Z^Z"" 


MUB 2 




X^ 


X^X'' 


MUB 3 


yA 


yB 


yAyB 


MUB 4 


X^Z"" 


yA^B 


^AyB 


MUB 5 


^AyB 


yA^B 


Z^X"" 



surements ll36l[37ll . Let us briefly review MUB, their proper- 
ties, and physical importance in the context of quantum mea- 
surement. 

Assume that {|ai)}^^o^ and {|&i)}f^Q^ are two different ba- 
sis sets for the d-dimensional Hilbert space 7Y. They are cafled 
mutually unbiased if they fulfill the following condition: 

As an example, for d = 2 (the case of a single qubit) it is 
easy to verify that the eigenvectors of the three Pauli matri- 
ces, X, Y, Z, denoted respectively by {|±)}x, and 
{|0), |1)}, constitute a set of pairwise MUB. In general, the 
maximum number of MUB for an arbitrary dimensional vec- 
tor space is not yet known, however, it has been proved that it 
cannot be greater than d + 1. In addition, for d being a power 
of prime, it has been proved that the number of MUB is ex- 
actly d + 1 and explicit construction algorithms are already 
known IjtIIssIi . For the case of n-qubit systems {d = 2") one 
can show the set of 4" — 1 Pauli operators, Ek = ®^i=iE]^(^^ ^-j , 
where Ea S X, Y, Z}, can be partitioned into 2" + 1 dis- 
tinct subsets, each consisting of 2" — 1 mutually commuting 
observables. All the operators in each subset have a set of joint 
eigenvectors. The eigenvectors of all subsets then form MUB 
ll39il40ll . Table Uillustrates such a MUB based partitioning of 
the 2-qubit Pauli operators. 

The importance of MUB which is relevant to our discus- 
sion is their application in quantum state estimation. To deter- 
mine the density matrix of a d-dimensional quantum system 

— 1 independent real parameters must be determined. The 
most informative (sub)ensemble measurements of an observ- 
able n of the system (whose spectrum is non-degenerate) pro- 
vide d — 1 independent data points, namely the probabilities 
Tr(p7ri), where ft = J^i'^i'^i is the spectral decomposition 
of 17 with spectrum uji.^ Thus, to fully determine the density 
matrix we must measure at least (d^ ~ 1) /{d — I) = d+1 dif- 
ferent noncommuting observables. In this sense measurement 
of the observables corresponding to MUB is optimal, because 
this requires the smallest possible number of noncommuting 
measurements. Moreover, due to the finiteness of ensembles 



^ In the non-degenerate case there are d — 1 orthogonal projectors ni since 
the Hilbert space is d-dimensional and Y2i ''^i = ^■ 



any repeated measurement will give rise to statistical errors. 
Naturally, to reduce such errors one must increase the size of 
ensembles and then repeat the measurements. However, it has 
been shown that a set of d + 1 MUB measurements provides 
the optimal estimation of an unknown quantum state, i.e., gen- 
erates minimal statistical error (if such MUB exist) iIjtIi . 

Now, we demonstrate that MUB measurements for state to- 
mography yields another version of AAPT. Let us first spe- 
cialize to the single qubit case. As noted earlier, to deter- 
mine a general quantum dynamical map on a single qubit us- 
ing AAPT, one attaches an ancilla and performs quantum state 
tomography at the end. In this case, the dimension of the com- 
bined Hilbert space is d = 2^ (we assume that the dimension 
of the ancilla is the same as that of the system). It follows 
from the general arguments above that one can use d + 1 ~ 5 
MUB measurements to determine the final state of the com- 
bined system (see Fig. |3]l. These measurements are in fact 
optimal in the sense explained earlier. The first (as always, 
ensemble) measurement provides four independent outcomes 
and each of the remaining ones yields three independent re- 
sults, which totals, as required, 4 + (4 x 3) = 16 results. 

It should be noted that even if the local state of the ancilla 
is known (i.e., if we know the expectation values of I^, 
I'^^X^,!^® Y^ and ® Z^ from prior knowledge about 
the preparation and trace-preserving property of the quan- 
tum map), the number of required measurements is still five 
ll4ll I42I1 . This can easily be seen from Table U For a non 
trace-preserving map we need five measurements of the (com- 
muting) operators of the first and the second columns (the ele- 
ments of the third column are products of the operators in the 
first two columns). If we know the local state of the ancilla 
B, the first three measurements of the second column are re- 
dundant. However, since the operators in the first column do 
not commute we still need to perform three (ensemble) mea- 
surements corresponding to the first three rows. The remain- 
ing two measurements related to the fourth and fifth rows are 
also necessary and they correspond to measuring the correla- 
tions of the principal qubit and the ancilla. Thus, the overall 
number of required MUB measurements in the case of trace- 
preserving maps is still 5. This argument is independent of the 
basis chosen, because in any other basis, due to noncommuta- 
tivity of the Pauli operators, the measurements corresponding 
to the local state of the ancilla always appear in different rows. 

For the case of n-qubit AAPT, the dimension of the joint 
system-ancilla Hilbert space is d = 2^". In this Hilbert space 
four different strategies can be devised: (i) using 16" (sep- 
arable) joint single-qubit measurements on the n-qubit sys- 
tem and the n-qubit ancilla (as explained earlier — Fig. |2|, 
(ii) using 5" MUB based measurements (tensor products of 
MUB based measurements of two-qubit systems), (iii) using 
d + 1 = 4" + 1 MUB based measurements on all 2n qubits, or 
(iv) using different combinations of single-, two-, and multi- 
qubit measurements including MUB based measurements (the 
number of measurements ranges from 4" + 1 to 16"). In what 
follows, we focus on method (iii) because it is the most eco- 
nomical in terms of the total number of measurements. The 
main drawback of performing a MUB based measurement on 
all 2n qubits is that it requires many-body interactions be- 
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FIG. 3: Schematic diagram of non-separable A APT. In this scheme 
the joint separable measurements of Fig. |2] have been replaced by 
(collective) mutually unbiased bases measurements on the two sys- 
tems. 



tween all 2n qubits. From an experimental point of view, 
such many-body interactions are not naturally available. This 
does not mean that they cannot be simulated, but as we will 
see this comes at a high resource cost. This is a strong re- 
striction which seriously affects the advantage of method (iii). 
According to our earlier discussion, the general multi-qubit 
observables in a MUB based measurement are generated from 
2^" + 1 noncommuting classes (or partitions) of 2ri-qubit op- 
erators {[iJi], . . . , [i?22"+i]}, where each class [Ek] contains 
2^71 _ 2 commuting observables, and Ek = ®'i=iEl^^^ with 
El^ e {I.X'-^, Z^} (for the case of three qubits refer to 
Fig. 2ofRef. H). 

In principle, one can simulate such many-body interactions 
from single- and two-qubit gates (e.g., CNOT). We next ar- 
gue that the complexity of such a quantum simulation scales 
at least as 0{n'^) or 0{'n?) depending, respectively, on the 
availability of non-local or local MUB measurements. 

All operators Em and £"„ that belong to the same class 
[Ek] commute and are composed of tensor products of iden- 
tity and/or Pauli operators. However, they cannot be simulta- 
neously measured locally, i.e. by using only single-qubit ob- 
servables. The reason is that each local measurement i?* , ^ 

for the operator Em completely destroys the outcome of mea- 
suring i?^!^^ for the other operator £"„, due to noncommuta- 
tivity of the PauU operators. It is simple to see that for the non- 
separable measurement of an operator such as . . . Z^", 
we need 2n sequential CNOT gates. To measure a more gen- 
eral observable such as Ek ^ E\^. . . -E^"^, where S^. G 
{J, X'\ Y^, Z^}, we need 0{n) additional single-qubit rota- 
tions to make appropriate basis changes. Therefore, for mea- 
suring n such general operators from the class [Ek], one needs 
to realize 0{n^) basic quantum operations. The condition 
for such a construction is the possibility of addressing arbi- 
trary distant pairs of qubits (i.e., having access to non-local 
two-body interactions). This is an important point, because 
in practical realizations the spatial arrangements of the qubits 
or other technological reasons may limit the interactions be- 
tween distant qubits. If only nearest neighbor gates can be im- 
plemented then pairs of qubits must be brought close to one 
another (e.g., via swap gates), which incurs a cost of 0{n) 



operations per pair [43]. In this case, we need 0{n'^) quan- 
tum gates to simulate the required multi-qubit measurements. 
It should also be noted that in such a simulation the scaling of 
execution time and possible (operational) errors in the mea- 
surements will introduce additional experimental complica- 
tions. 



C. Generalized measurement 

In principle, it is also possible to perform the required quan- 
tum state tomography at the output of AAPT by utilizing only 
a single generalized measurement or POVM ll44l l45ll . Sup- 
pose that we want to determine an unknown state p of our 
quantum system. In a d-dimensional Hilbert space charac- 
terization of p requires determination of cP — 1 independent 
real parameters. In order to design a scheme for determina- 
tion of /9 by a single quantum observable, we need to attach a 
(i' -dimensional ancilla (B) with a known initial state r to our 
principal system (A). In the scheme proposed in Ref. liill . 
one should measure one of the observables of the combined 
system (AB), a "universal quantum observable". 



dd' 

a=l 



Pn 



(11) 



where is a normal operator and the spectrum Aq should 
be non-degenerate such that the projections Pa constitute a 
complete set of dd! — 1 commuting observables. Since the 
projections all commute, one can measure all of them si- 
multaneously using a single apparatus. Such (repeated en- 
semble) measurements provide us with dd' — 1 probabilities 
Pa = Tr(Pap ® r).^ The dimension of the ancilla must be 
greater than or equal to the dimension of the system, d' ^ d. 
If we take p = j^nm Pnm\n) {m[ and r = J2a[3'^ai3\a){(3[, 
then we obtain the following linear relation: 



(12) 



where -A/^^^ = 'rai3{mP[Pa[na). When d = d' and the 
measured observable couples A and B in a manner such that 
Ma^mn = ^'^mji invcrtiblc (detAf 7^ 0), a linear inversion 
can reveal the unknown state p |47]. 
To be specific, we choose as follows: 



n = J2aE^ 

a=l 



•E^ 



(13) 



where ({^-^jf^i) is a set of orthonormal basis op- 

erators for the space of linear operators on Ha {'Hb)-'' Using 
the representation T = J2a Tr(I'£'J)-Ba (for any operator T), 
it is not hard to see that the ensemble average of an arbitrary 
operator O (on Ha) is equivalent to an ensemble average of 
the following function of fl: 



F - V Tr(0£;„^t) 



E' 



* In fact, as noted in l47ll . the set of operators {Trs(rPa)} constitutes in 
TiA a minimal informationally complete POVM |'46:]. 

^ The operators {Ea } should be normal: [Ea, eI] = 0, which makes them 
observable in the sense defined in Ref. l44ll . In the multi-qubit case the 
basis operators can be taken as tensor products of the Pauli operators. 
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FIG. 4: Schematic diagram of a QPT by using POVM. Here we have 
used the idea of "universal quantum observable" (44ll . To accom- 
plish complete process tomography, one needs two more ancillas Bi 
and B2 (in addition to the one used in AAPT, B) and two universal 
quantum observable Qi and ^2. 



on p (g) r, i.e., 

{0)p = {Fomp^r- (14) 

Therefore estimation of the ensemble average (0)p of an op- 
erator O acting on the principal system A, can be achieved by 
measuring Fo{^) on the joint A and B system. This allows 
for the estimation of every ensemble average for the principal 
quantum system. 

The above general scheme can also be utilized for quantum 
process tomography (Fig. 3 in Ref. liill ). It is sufficient to 
consider the AAPT scheme and attach two additional ancillas 
(one for the system and another for the ancilla of the AAPT 
scheme), and then measure jointly two universal observables 
(Fig. nil. In this manner, to characterize the dynamics of n 
qubits, the number of required ancillary qubits increases from 
n (in AAPT) to at least 3??. This can be easily understood 
via a simple counting argument. In order to extract complete 
information about a quantum dynamical map on n qubits (en- 
coded by 2''" independent parameters of the superoperator x) 
in a single measurement, one needs a Hilbert space of dimen- 
sion at least 2''" on which the information can be imprinted 
unambiguously. 

There are two major disadvantages in using such a POVM 
compared to all other QPT schemes. First, the POVM scheme 
requires a general many-body interaction between all 2n 
qubits that are measured through each 17^. This interaction 
cannot be efficiently simulated, i.e., it requires an exponential 
number of elementary single- and two-qubit quantum gates. 
Indeed, the above universal quantum observable scheme is 
very difficult to implement in practice, because it implies mea- 
suring an observable il (or a commuting set {Pa}) which thor- 
oughly entangles the system and the ancilla(s). There is an 
alternative method to implement the above scheme ll47ll . by 
interacting system and ancilla for a specific time duration t 
through a known unitary operator U (or known Hamiltonian 
H), and then measuring the simplest possible non-degenerate 
observable fi, namely a factorized quantity = lu^ (g) cu^ 
ll44i I47I1 . However, even this method still requires a many- 
body interaction (through U) which is difficult to prepare. The 
operator fl has maximal Schmidt number and generally can- 
not be simulated using a polynomial number of elementary 
gates. It is known that, in general, 0(4^) elementary single- 
and two-qubit gates are necessary to simulate many-body op- 



erations acting on A'^ qubits ll48ll (see also Ref. jsill for differ- 
ent measures of complexity of a given quantum dynamics, and 
Ref. 1I49I1 for the concept of entanglement cost of a POVM). 

V. DIRECT CHARACTERIZATION OF QUANTUM 
DYNAMICS 

Recently a new scheme for quantum process identification 
was proposed and termed "direct characterization of quan- 
tum dynamics" (DCQD) JItI H. It differs in a number 
of essential aspects from SQPT and AAPT. In DCQD, sim- 
ilarly to AAPT, the degrees of freedom of an ancilla system 
are utilized, but in contrast it does not require inversion of 
a full (P X (P matrix (hence "direct"); it requires different 
input states; and uses a fixed measurement apparatus (Bell 
state analyzer) at the output. The main idea in DCQD is to 
use certain entangled states as inputs and to perform a sim- 
ple error-detecting measurement on the joint system-ancilla 
Hilbert space. A combination of these input states and mea- 
surements give rise to a direct encoding of the elements of the 
quantum map into the measurement results, which removes 
the need for state tomography. More precisely, by "direct" 
we mean that the measured probability distributions (on an 
ensemble of the setting) are rather directly related, i.e., with- 
out the need for a complete inversion, to the elements of x- 
In essence, in DCQD the % matrix elements of linear quan- 
tum maps become directly experimentally observable. For the 
case of single qubit, the measurement scheme turns out to be 
equivalent to a Bell-state measurement (BSM). In DCQD the 
choice of input states is dictated by whether diagonal (popula- 
tion) or off-diagonal (coherence) elements of the superopera- 
tor are to be determined. Population characterization requires 
maximally entangled input states, while coherence character- 
ization requires non-maximally entangled input states. In the 
following, we review the scheme for the case of qubits. For a 
generalization of the scheme to higher-dimensional quantum 
systems see Ref. ifisll . 

Let us consider the case of a single qubit and demonstrate 
how to determine all diagonal elements of the superoperator, 
{Xmm}, in a single (ensemble) measurement. We choose 
{I , , Z"^} as our error operator basis acting on the 

principal qubit A. We first maximally entangle the two qubits 
A (the principal system) and B (the ancilla) in a Bell-state 
|$+) = (|00) + |ll))yiB/\/2(aninstance of a stabilizer code), 
and then subject only qubit A to the map £. 

A stabilizer code is a subspace He of the Hilbert space of 
71 qubits that is an eigenspace of a given Abelian subgroup S 
(the stabilizer group) of the 71-qubit Pauli group, with eigen- 
value +1 lUllsO]. In other words, for every |^'c) G He and 
all Si e S, we have S'^l^'c) = |*c>, where the Si's are the 
stabilizer generators and [Si, Sj] ~ for all i and j. Consider 
the action of an arbitrary error operator E on the stabilizer 
code state: E\'^c)- The detection of such an error will be 
possible if the error operator anticommutes with (at least one 
of) the stabilizer generators: {Si,E} = 0. To see this note 
that 

Si{E\^c)) = ~E{Si\^c)) - -(£^l*c)), 
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i.e., E\^c) is a —1 eigenstate of Si. Hence measurement of 
Si detects the occurrence of an error or no error (—1 or +1 
outcomes, respectively). Measuring all the generators of the 
stabilizer then yields a list of errors ("syndrome"), which al- 
lows one to determine the nature of the errors unambiguously. 

The state |$+) is a +1 eigenstate of the commuting oper- 
ators Z^Z^ and X^X^, i.e., it is stabilized under the ac- 
tion of these stabilizer generators. It is easy to see that any 
non-trivial error operator Ei G {I,X^, Y^, Z"^} acting on 
the state of the qubit A anticommutes with at least one of the 
stabilizer generators, and therefore by measuring them simul- 
taneously we can detect the error: 



A vB 



X^'X 

z^z^ 



{E^\^+))^±{Ef\^+)). 



Note that measuring the observables Z Z and X X is 
indeed equivalent to a BSM, and can be represented by the 
four projection operators P^± = I^I**) and P,j,± = 

where = (|00) ± |ll))/\/2 , and = 

(|01) ± |10))/\/2arethe Bell states. The probabilities of ob- 
taining the no-error outcome /, bit-flip error X"^, phase-flip 
error Z^, and both phase-flip and bit-flip errors Y^, on qubit 
A can be expressed as: 



(15) 



where m = 0,1,2,3, and the projectors P,„, for ra = 
0, 1, 2, 3, correspond to the states $+, v]/-, and re- 
spectively. Here £{p) is a shorthand for {£ ® I){p)- Equa- 
tion ( fTSl l is a remarkable result: it shows that the diagonal 
elements of the superoperator are directly obtainable from an 
ensemble BSM. This is the core observation that leads to the 
DCQD scheme. In particular, we can determine the quantum 
dynamical populations, Xmm, in a single ensemble measure- 
ment (i.e., by simultaneously measuring the operators Z^Z^ 
and X^X^) on multiple copies of the state |$^)). 

To determine the coherence elements, Xm^n, a modified 
strategy is needed. As the input state we take a non-maximally 
entangled state: |$+) a|00) + with |a|, ^ 

{0, l/\/2} and Im(a/3) ^ 0. The sole stabilizer of this state 
is Z^Z^ . The spectral decomposition of this stabilizer is 
Z^Z^ — P_i_]^ — P_i, where P±i are projection operators 
defined as P^i = + P^- and P_i = + P<^- . Now, 
it is easy to see that by measuring Z^Z^ on the output state 
£{p), with p = we obtain: 



and 



Ti[P+i£{p)] = xoo + X33 + 2Re{xo3){Z^), (16) 



Tt[P_i£{p)] = xii + X22 + 2lm{xi2){Z^), (17) 



where (Z^) = Tr{pZ^) ^ because of our choice of a non- 
maximally entangled input state (|a|, |/3| ^ {0, l/^/2}). The 
experimental data, 1r:[P±i£{p)\, are exactly the probabilities 
of no bit-flip error and a bit-flip error on qubit A, respec- 
tively. Since we already know the Xmm's from the popula- 
tion measurement, we can determine Rc(xo3) and Im(xi2). 



A- 

B- 



BSM 



FIG. 5: Schematic diagram of the DCQD scheme. The system and 
the ancilla are prepared in one of the input states as in Table Ull and 
after subjecting the system to the map £, the combined system is 
measured in the Bell-state basis. 



After measuring Z^Z^ the system is in either of the states 
p±i = P±i£{p)P±i/Ti[P±i£{p)]. Next we measure the 
expectation value of a normalizer operator A'^, for example 
X^X^, which commutes with the stabiHzer Z^Z^.^ We 
then obtain the measurement results 

Ti[Np+,] = [{xoo-X33){N)+2tlm{xo3){Z^N)]/Tr[P+i£{p)] 



and 



Tr[7Vp_i] = [ixii-X22){N)~2iReixi2){Z^N)]/Ti[P^,£{p)], 

where {Z^), (N), and (Z'^N) are all non-zero and already 
known. In this manner, via a simple linear algebraic calcu- 
lation, we can extract the four independent real parameters 
needed to calculate the coherence elements xo3 and xi2- It 
is easy to verify that a simultaneous measurement of the sta- 
bilizer, Z^Z^, and the normalizer, X"^X^, is again nothing 
but a BSM. However, in order to construct the relevant infor- 
mation about the dynamical coherence, we need to calculate 
the expectation values of the Hermitian operators P$+ ± P$~ 
and P^+ ± P^-. 

In order to acquire complete information about the coher- 
ence elements of the unknown dynamical map £, we perform 
an appropriate change of basis by preparing the input states 
H^H^\^+) and S^S^ H^H'^\^+), which ai-e the eigen- 
vectors of the stabilizer operators X^X^ and Y^Y^ . Here 
H and S are single-qubit Hadamard and phase gates acting on 
the systems A and B. At the output, we measure the stabilizers 
and a corresponding normalizer, e.g., Z^Z^ , which are again 
equivalent to a standard BSM, and can be expressed by mea- 
suring the Hermitian operators P$+ ± P^+ and P$- ± P^^ 
(for the input state H^H^\^X))^ andP$+ iPvp- andP$-± 
(for the input state S^S^ H^H^\^X))- Figure H illus- 
trates the DCQD scheme. 

Overall, in DCQD we only need a single fixed measurement 
apparatus capable of performing a Bell state measurement, for 
a complete characterization of the dynamics. This measure- 
ment apparatus is used in four ensemble measurements each 
corresponding to a different input state. Figure |5] and Table Ull 
summarize the preparations required for DCQD in the single 
qubit case. This table implies that the required resources in 



A normalizer operator A'^ is a unitary operator that preserves the stabilizer 
subspace but is not in S. The normalizer group M commutes with the 
stabilizer group <S. 
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TABLE II: One possible set of input states and measurements for di- 
rect characterization of quantum dynamics (x) for a single qubit. 
Here !$+) = a|00) + {\a\ / 0, l/x/2), \<!>i)x(Y) = 

"1 + +)x(Y) -)x(Y) (kl / 0, 1/V2 and Im(a^) / 0) 

and {|0), {i±)jsr}, {i±)y} are eigenstates of the Pauli opera- 
tors Z, X, and Y. The fourth column shows the BSM measurement 
equivalent to stabilizer -I- normalizer measurements. 



input state 



1* + ) 



Stabilizer 



Measurement 
Normalizer BSM 



Z^Z'' 

yAyB 



N/A 



± P P 



± P, 



± P,T, 



± P P 



± P, 



output mn 



00,11.22,33 
03,12 
01,23 
02,13 



DCQD are as follows: (a) preparation of a maximally entan- 
gled state (for population characterization), (b) preparation of 
three other (non-maximally) entangled states (for coherence 
characterization), and (c) a fixed Bell-state analyzer. 

Our presentation of the DCQD algorithm assumes ideal 
(i.e., error-free) quantum state preparation, measurement, and 
ancilla channels. However, these assumptions can all be re- 
laxed in certain situations, in particular when the imperfec- 
tions are already known. A discussion of these issues is be- 
yond the scope of this work and will be the subject of a future 
publication liSLI . 



VI. DISCUSSION AND RESOURCE COMPARISON 

In this section we present a discussion and comparison of 
the various QPT schemes described in the previous sections, 
and highlight the important features of each scheme, as illus- 
trated in Tables Ullland llVI The goal is to provide a (physical) 
resource analysis and guide for choosing the appropriate QPT 
scheme, when available resources and the particular system of 
interest are taken into consideration. 



A. Scaling of the Required Number of Experimental 
Configurations witli tlie Number of Qubits 

For characterizing a quantum dynamical map on n qubits 
we usually perform measurements coiTesponding to a tensor 
product of the measurements on the individual qubits. An 
important example is a quantum information processing unit 
with n qubits. DCQD requires a total of 4" experimental con- 
figurations for a complete characterization of the dynamics, 
where the total number of experimental configurations is de- 
fined as the number of input states times the number of non- 
commuting measurements per input — see Table |III] This is 
a quadratic advantage over SQPT and separable AAPT, both 
of which require a total of 16" experimental configurations. 
However for quantum systems without controllable two-qubit 
operations, implementation of the DCQD scheme is hard, here 
SQPT is the most efficient scheme. 

In principle, the required state tomography in AAPT could 



also be realized by non-separable (global) quantum measure- 
ments. These measurements can be performed either in the 
same Hilbert space, with 4" + 1 MUB measurements, or in 
a larger Hilbert space, with a single generalized measure- 
ment. Both methods require many-body interactions among 
2?? qubits, which are not naturally available. For the AAPT 
scheme using MUB measurements, one can simulate the re- 
quired many-body interactions using a quantum circuit com- 
prising 0{n^) [0{n^)] single and two-qubit quantum elemen- 
tary gates, under the assumption of available non-local [lo- 
cal] two-body interactions. On the other hand, in DCQD the 
only required operations are Bell-state measurements, each of 
which requires one CNOT and a Hadamard gate. This results 
in a linear, 0{n), scaling of necessary quantum operations for 
realization of each experimental configuration in DCQD — see 
TableHV] 



In general, in the 2^"-dimensional Hilbert space of the 2n 
qubits of the system and the ancilla, one could devise in- 
termediate strategies for AAPT using different combinations 
of single-, two-, and many-body measurements. The num- 
ber of experimental configurations in such methods ranges 
from 4" + 1 to 16", which is always larger than that of 
DCQD, which requires 4" BSM setups. Therefore, in the 
given Hilbert space of n qubits and n ancillas, DCQD requires 
fewer experimental configurations than all other known QPT 
schemes. In this sense, DCQD has an advantage over AAPT 
in a Hilbert space of the same dimension. Moreover, using 
DCQD one can transfer logj 2^" bits of classical information 
between two parties, Alice and Bob UtIi . which is optimal 
according to the Holevo bound |[lt]. This is a similar context 
to the quantum dense coding protocol |[l|]. Alice can realize 
the task of sending classical information to Bob by applying 
one of 2^" unitary operator basis elements to the n qubits in 
her possession and then send them to Bob. Bob can decode 
the message by a single measurement on his 2n qubits using 
the DCQD scheme. In other words, the total number of possi- 
ble independent outcomes in each measurement in DCQD is 
2^", which is exactly equal to the number of independent de- 
grees of freedom for a 2?i-qubit system. Therefore, a maximal 
amount of information can be extracted in each measurement 
in DCQD, which cannot be improved upon by any other pos- 
sible QPT strategy in the same Hilbert space. 



For characterizing the dynamics of n qubits in a single 
generalized (POVM) measurement unambiguously, a Hilbert 
space of dimension at least 2^" is required. In order to im- 
plement such a POVM, one needs to realize a global normal 
operator (a single universal quantum observable) acting on the 
joint system-ancilla Hilbert space, of the form of Eqs. (fTTI) 
and ( fT3] l. Such generic operators cannot be simulated in a 
polynomial number of steps. It is known fll] that in gen- 
eral at least 0(4^") single- and two-qubit basic operations are 
needed to simulate such general many-body operations acting 
on 2n qubits. 
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TABLE III: Required physical resources for the QPT schemes: Standard Quantum Process Tomography (SQPT), Ancilla-Assisted Process 
Tomography (AAPT) using joint separable measurements (JSM), using mutual unbiased bases measurements (MUB), using generalized mea- 
surements (POVM), and Direct Characterization of Quantum Dynamics (DCQD). 



Scheme 



dimCH)" 



Nin, 



as ./input 



measurements 



required interactions 



SQPT 
AAPT 

DCQD 



JSM 

IVIUB 

POVM 



2" 
22„ 



4" 
1 
1 
1 

4" 



4" 
16" 
4" + 1 
1 
1 



16" 
16" 
4" + 1 
1 

4" 



1-qubit 
joint 1-qubit 
MUB 
POVM 
BSM 



"Ti.: the Hilbert space of each experimental configuration 
''total number of noncommuting measurements per input 
"^total number of experimental configurations = 



inputs 



X N, 



single-body 
single-body 
many-body 
many-body 
single- and two- body 



meas7input 



B. Accuracy Considerations 

Due to the finiteness of the number of measurements that 
can be performed in practice when estimating an ensemble 
average, it is evident that estimation of an unknown quantum 
map through any of the QPT schemes gives rise to some er- 
ror. Such statistical errors can in principle be reduced by in- 
creasing the size of ensembles. A relevant question in QPT 
discussions is then how the accuracy of estimations in differ- 
ent QPT schemes depends on the ensemble size (N). Finite 
size scaling behavior of this accuracy (or error) can provide 
another practical figure-of-merit for comparison of different 
QPT schemes. Here, our discussion is just tangential and very 
incomplete so that it just aims at showing just a rough picture 
of the issue. A complete investigation of the finite-size errors 
is not our goal in this paper. Another issue that we partially 
address here is numerical error due to the inversion required 
in some QPT schemes. 



1. Finite Ensemble-Size Effects 



There is a huge literature regarding analysis 
timation errors or quantum statistics 11521 l53l \54 



uantum es- 




58 , m MMMMMM 

7^ "vmU%mV^ Our aim here is to give 

a very brief discussion of estimation errors in different QPT 
schemes through a special example. At the end of this sub- 
section we go a bit further and provide a sketchy discussion 
of more standard figures-of-merit. However, a more complete 
investigation of this subject is beyond the goals of this paper 
and needs a separate study per se. 

In all QPT schemes measurements are performed of one or 
more observables {O^^'j^j^, where X denotes the scheme: 
X € {SQPT, AAPT-JSM, AAPT-MUB, AAPT-POVM, 
DCQD}. E.g., /^AAPT-JSM = 16" (all operator basis for the 
entire Hilbert space of system and ancilla), i^AAPT-povM = 1 
(the universal observable fl), A'dcqd = 4" (4 Bell state mea- 
surements per principal qubit, one for the superoperator pop- 
ulation, three for the coherences — ^here it makes no difference 
if measurements commute). For notational simplicity let us 
omit the (X) superscript. Each observable (given scheme X) 

has a spectral decomposition: Ok — J2i=o^ aI'^'p/*^', where 



(k) (k) 

are the eigenvalues and are projection operators. 
The number i/^ of distinct projectors is the number of possi- 
ble measurement outcomes for a given observable Ok (which 
is typically the dimension of the relevant Hilbert space). E.g., 
in AAPT-POVM (where there is only a single observable), 
= 16", and in DCQD i/fe = 4" for all k (n-fold tensor 
product of Bell state measurements on qubit pairs). We can 
also interpret i/k as the dimension of the probability space 
associated with a random variable Yk that can take values 

ie{0,...,i^k-i}- 

Given an observable Ok, we must be able to unambiguously 
determine the index i of which projection operator (or eigen- 
value) was measured. For example, when we measure the 
Pauli operator Z, the projectors (eigenvalues) are Pq = |0) (0| 
(Aq = 1) and Pi = |1)(1| (Ao = —1), and we must have 
a device (e.g., a Stern-Gerlach detector) which unambigu- 
ously reveals whether the final state has spin up (Pq) or down 
(Pi). In other words, the raw experimental outcomes are de- 
tector clicks in bins that count how many times nf''' each in- 
dex i has been obtained. The resulting empirical frequencies 



{/^''^ EE n\'^'/Nk}, where Nk I]i=o " approx- 
imations to the true probabilities {p[''^} of detector clicks: 

(fc) _ Tr[£{p)P^''^]. In terms of the random variable de- 

[k) 



(k) 



-1 (fe) 



P: 



scription mentioned above, we have Pr(Yfc = i) = p. 

For a given observable Ok, repetition of the experiment or 
increase in the sample size Nk can reduce the error in the 
probability inference A^"* = \p^^ — |. However, we 
are in general interested in the expectation values of the ob- 
servables Ok, that can be obtained from the probability dis- 
tribution {pf }£o^- I-e- 

we would like to know the true 

mean (Ok) = Tr[£(p)0,] = ^^=0 ' 'pf \ which we 
estimate using the empirical frequencies to get the empiri- 
cal mean fik = J2ito^ '^i''^ fi'^^ ■ The central Umit theo- 
rem Issll (or the Chernoff inequality jsill ) states that in the 
limit Nk oo the probability that the empirical mean Hk 
is far from the expected value (Ok), is very small. More 
precisely, defining the true standard deviation as usual as 

'^fe ^ y(oiy~w [where (o^) ^ j::^\x^^)^p^\ 

and letting 2^/2 be the cutoff for the upper tail of the nor- 
mal distribution N((Ofc), cr^) having probability a/2, we have 
asymptotically: Pr(|(Ofe) — fik\ ^ ~T0°) = 1 — «• Here a 
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TABLE IV: Resource analysis of the QPT scheme (for the case in which the probabilities {p,-''' } are distributed uniformly): Standard Quantum 
Process Tomography (SQPT), Ancilla- Assisted Process Tomography using joint separable measurements (JSM), using mutual unbiased bases 
measurements (MUB), using generalized measurements (POVM), and Direct Characterization of Quantum Dynamics (DCQD). 

Scheme A^c^p." 1-qubit gates/meas. 2-qubit gates/meas. -/Vgaics/mcas. -^ovcnUl^ 

SQPT 16" 0{n) N/A 0(n) 0{nl6'^) 

JSM 16" 0(n) 0(n) 0(n) C'(nl6") 

AAPT f^fjg 4" + l 0(ri=) O(n^) [©(n^)] Oin') [Oin'^)] 0(n^4") [0(n='4")] 

POVM 1 0(4^") 0(4^") 0(4^") 0(4^") 

DCQD 4" 0{n) 0{n) 0(n) 0{n4") 

"as defined in Table lllll 

'overall complexity = A^exp. X A'gates/meas. 



represents the confidence interval. This result allows us to 
compare the estimates of any two means, by equating their 
confidence intervals. By replacing the tr ue standard devia- 
tion by the empirical one, i.e., by = y/ {fJ''^)k — where 

{fi'^)k = X]r=o ^('^i'^^)^//''^ (an excellent approximation in 
large samples), we have 



AX') 



(18) 



as the criterion for having a confidence interval of equal length 



around the two sample means fi)^ 



and fi^i^ 



I.e., to contain 

the unknown true means (Ok) and {Ok') with equal proba- 
bility. Here we have reintroduced the QPT label (superscript 
X) to stress that this criterion holds for the comparison of es- 
timates of any two means, across both k and X. This result 
shows that, assuming the standard deviations do not scale with 
i^fe, equal confidence in estimates of two expectation values of 

simply requires equal sam- 



two observables O). and O)., 



pie sizes iV^^^ and A^^f 

However, let us note that the above statistical argument is 
rigorous only in the limit Nk oo. That the situation is dif- 
ferent for finite sample sizes can be appreciated via the follow- 
ing examples, for which we first recall the Chernoff inequality 
ll84ll . The version of this inequality which is best suited to our 
present purpose is as follows. Assume that an event 7 occurs 
with the true probability ^(7). We estimate this probability 
by performing N independent trials. The inferred probability 
is then pn{i) = nM{'^)/N, where nM{'^) is the number of 
occurrences of 7 in the trials. Then for any A G [0, 1] the 
Chernoff inequality is 

Pr {\pn{i) - P{l)\ ^ ^P{l)) ^ e-P(^)^'^V3. (19) 



An immediate result of this inequality is the following. Let 

(20) 



^(7;A,e) = ^4T7log-• 



p(7)A2 

Then for any A, e G [0, 1], if iV ^ iV(7; A, e), then with 
probability greater than 1 — e we have \pn{i)—p{i)\/p{i) ^ 
A. Roughly, if we wish Pn{i) to be within an error of at 
most A from ^(7), this can happen with a probability greater 



than 1 — e (for some e) when we perform at least A^(7; A, e) 
trials. It follows that a highly accurate estimation (A, e — ^ 0) 
requires many (TV — > 00) trials. In standard statistical error 
analysis A is usually taken to be the standard deviation cr at or 
at most 2(7 AT. 

From Eq. ( |20l i it is evident that if the probabilities {pi\ does 
not depend on the dimension of the Hilbert space (n) the num- 
ber of repetitions to fulfill an error e would not either — this 
number will only have a logarithmic dependence on the er- 
ror e. This implies that there are cases in which the statistics 
can be built up with a constant overhead in ensemble size- 
up to the logarithmic dependence on the error This can have 
highly useful and efficient applications for QPT in such cases. 
Nonetheless, it would be important to point out an intricate 
pitfall in (incautious) too general conclusions. To this aim, 
here we want to analyze a somewhat pathological example in 
which efficiency cannot be concluded from Eq. ( l20l i. Let us 
assume that we are dealing with fairly uniform probability dis- 
tributions {pP''}J'4o^ and compare two situations: tossing a 
coin (two possible outcomes: vi ~ 2) and estimating the prob- 
ability distribution of a random variable with V2 = 10^" dif- 
ferent possible outcomes. In the case of the coin it is clear that 
after N = 10^ tosses we will have a pretty good idea about the 
probabilities andp^"* ofheads vs tails. On the other hand, 
for the other random variable, after = 10^ measurements 
we will have not yet sampled the entire space of possible out- 
comes, so will not have been able to gather statistics represen- 
tative of all probabilities p] (some outcomes will not have 
ever occurred, thus their probabilities cannot be estimated). 
Consequently, we will not be able to accurately estimate any 
means. However, from Eq. ( l20l i. with the uniform probability 
distribution assumption p^'^-' = \/vk, we obtain 



Nk ^ Af(A,e) 



3^ log; 



i^feC(A,e) 



(21) 



In other words, for accurate estimation of means in the case of 
fairly uniform probability distributions it is sufficient to have 
Nk 5^ C(A, €)vk, where the prefactor C(A, e) encompasses 
both the estimation error A and the probability 1 — e to achieve 
that error We call condition (I2TI 1 the "good statistics" con- 
dition. It is important to note that this conclusion depends 
strongly on the assumption of fairly uniform probability distri- 
butions. Indeed, consider the case where the random variable 
with 10^" different possible outcomes is very strongly peaked 
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at two values {ii, 42}- In this case it behaves effectively like a 
coin, and we do not need N ^ 10^°. 

We thus see that a comparison of the different QPT meth- 
ods on the basis of fixed mean-estimation error will depend 
strongly on the properties of the underlying probability distri- 
butions {{p\''^ }ito^}k=i- ^ thorough study of the properties 
of these probability distributions as a function of QPT method 
X is beyond the scope of this paper However, let us specu- 
late on what would happen if the assumption of fairly uniform 
distributions were to hold for all k and X. The total number 
of ensemble measurements becomes 



Kx 

fe=i 



inputs, A; ' 



(22) 



where N- 



'X) 
inputs, 



is the number of initial states needed per ob- 

servable k for a given QPT scheme {X), and Nj, is found 
from the good statistics condition (I2TI 1. with i^k replaced by 

vj. . We can read off the values of N^Jts k ^^'^ from 
the second and third columns of Table [ffl] respectively. The 
values of are as follows. In the case of X=SQPT and 
AAPT-JSM the observables are all one-dimensional projec- 
tors so that = 1 Vfc. In the case of X=AAPT-MUB there 
are 2n qubits (i.e., a 4"-dimensional Hilbert space) and one 
performs quantum state tomography at the output by measur- 
ing a set of noncommuting 4" + 1 observables of the MUB 
basis, where each member of the MUB basis has a spectral 
resolution over Vk = 4" — 1 independent projective measure- 
ments.^ We already noted above that = 4" and 16" 
for X=DCQD, and AAPT-POVM, respectively. We observe 
that, for fairly uniform distributions, n'^^ grows exponen- 
tially with respect to the number of qubits n for non-separable 
process tomography schemes, with AAPT-MUB and AAPT- 
POVM at a distant disadvantage due to the inherent depth of 
their quantum circuits for simulating many-body interactions 
in each measurement (see the fourth column of Table llVl i. 
Collecting the results above, however, it follows from Eq. (|22] | 
that the total number of ensemble measurements N^-^^ scales 
as 16" in all QPT methods, to within a factor C(^)(A, e). 



How would the number of ensemble measurements, N, 



'X) 



change if the distributions are sharply peaked? For separable 
schemes, e.g., SQPT and AAPT-JSM, this would not result 
in any difference, since we already have z^^ = 1. How- 
ever, for non-separable schemes this would lead to substantial 
reduction of measurements since we would be dealing with 
effectively fixed-dimensional probability distribution spaces, 
e.g., v^. = const., instead of an exponential function of the 
number of qubits. Hence the question of the properties of the 

probability distributions {{Pi'^''}i'=o indeed important 

and will be the subject of a future study. 



' One of the 4" + 1 observables has 4" independent outcomes, which gives 
(4" + 1)(4" — 1) + 1 = 16" outcomes, which is sufficient to fully 
characterize the superoperator. 



2. Discussion of Figure-of-Merit 

One of the standard approaches in quantum estimation 
and quantum statistics to address estimation errors is via 
the Cramer-Rao bound (CRB) iH IHIll, I?!- Following 
Ref. ifssll . the CRB can be described as follows. Let us as- 
sume that {X)rm} G are the true (real-valued indepen- 
dent) parameters of x that are supposed to be estimated from 
a measurement data set Vx obtained through the scheme X - 
we remove the superscript R in the sequel without any risk of 
confusion. The true negative logarithmic likelihood function 
of the system generating that true data is defined by 



Kx i^k-l 

EE' 

k=l i=0 



logp,' 



(k) 



(k) 

where ' is the number of times the outcome i is obtained 
from £k measurements of Ok (total of J^k measurements) 
and (n,|'^'') = v'-f' ^k (where ( ) is the quantum average). If 
X e K.^ is an unbiased estimate of i-^-, (x) = X, the 
covariance of the estimate cov(x) = ((x — — x)^) 

satisfies the following matrix inequality: 




^ 0, 



(23) 



where T is the Fisher information matrix defined as 



mn.m'n' 



ix) = {V^,^,logC^^^\ 



_ dlogpf"^ 9 log p. 
" ^~d^ By , , 

Provided that J-{x) is positive and invertible, Eq. ( l23l l gives 
the following well-known form of the CRB; 



(k) 



cov(x) ^ J" ^ix) 



(24) 



Taking the trace of both sides and noting that var(x) = 
Tr[cov(x)], one can also find a scalar form for this equation. 
Equation (l24l i means that for any unbiased estimator the er- 
ror is lower-bounded by the inverse of the Fisher information. 
The Fisher information matrix is indeed a measure of infor- 
mation about X that exists in the data Vx ■ The special feature 
and indeed the power of this bound is that it is independent of 
how the estimate is obtained, for T is independent of the esti- 
mation mechanism. For the case of single-parameter estima- 
tion, the CRB can always be achieved a symp totically by using 
maximum likelihood (ML) estimation l54l l74l [83ll . That is, 
as the amount of data increases the ML estimate approaches 
the true answer with the error bars equal to those given by 
the CRB. However, for multi-parameter estimation there is, in 
general, no optimal estimator that can achieve this bound. See 
Ref. ll74ll and references therein for more information about 
the CRB, its quantum version, and its application to quantum 
state estimation. The above discussion may suggest that the 
(inverse of the) Fisher information matrix can be taken as a 
good figure-of-merit for a quantum estimation process. How- 
ever, the very nature of independence from estimation method 
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means that the Fisher information matrix is not so useful for 
the purpose of comparing different QPT schemes — our goal 
in this paper. 

A more promising and physically-motivated approach, that 
justifies using the Chernoff bound for the purpose of quan- 
tum state/process estimation as we did earlier, has been pro- 
posed very recently, and is called the quantum Chernoff bound 
(QCB) U9, 80, 81, 82]. The physical interpretation of this 
quantity is as follows: assuming that we have access to all 
types of measurements — whether local or collective — on all 
ensembles, the QCB measures the error in distinguishing a 
state p from another state p. The probability of a wrong 
inference, i.e., mistaking p for p, has the asymptotic form 
Pe e^''^^™(''''^), where 

Acb(p,p)= min Tr[pV'""], (25) 

is the QCB , and ^ Acb ^ 1- The maximum is attained for 
P = P 11791 [sil [82I1 . Recently, the QCB has been considered 
as a natural figure-of-merit in evaluating the performance of 
different measurement scenarios for qubit tomography [81]. It 
also has been used for quantum hypothesis testing and distin- 
guishability between density matrices Isol [82ll . Considering 
the fact that a generic x matrix is formally in the category of 
density matrices, the application of the QCB can in principle 
be extended to QPT. That is, one can in principle calculate 

{x,x) for estimation of a quantum process x through 
any QPT scheme X and then take an average over all possible 
processes with a suitable probabiUty measure dfi{x) ll8ll[82ll . 
The average QCB 

a(^) = J dp{x)E^p{x\x)A^CB\x,x), (26) 

where p{x\x) is the probabiUty of estimating x given the true 
process x, may prove a more useful figure-of-merit. A more 
complete analysis, along with possible numerics, that explic- 
itly shows the performance of different QPT schemes (similar 
to the analysis of Ref. ifsill ). is yet to be performed. One im- 
portant point, however, is the issue that may be caused by the 
assumption of availability of all types of measurements (in- 
cluding collective measurements) in this bound and whether 
they are important in achieving the bound or not. This may 
in turn complicate usage of this tool as a completely suit- 
able figure-of-merit for a comparative study of different QPT 
schemes. For completeness, let us mention that a different 
investigation of physically good figures-of-merit (or distance 
measures) for quantum operations has also been performed in 
Ref. iH. 

Other characteristics of the QPT schemes may also play 
significant roles in the propagation of errors in the inferred 
quantum map £. Indeed, the effect of preparation, i.e., how 
different input states can affect efficiency of the estimation of 
unknown maps, must be explored as well - for a recent study 
see Ref. llssll . For the case of AAPT, as explained earlier, it 
is already known that using maximally entangled input states 
is favored, because they result in smaller experimental errors 
than separable states. For DCQD an analysis of how different 
input states affect performance of the estimation is underway 



ll23ll . Without a full understanding of the role of preparation, 
the scaleup of physical resources in different QPT strategies 
for finite ensemble sizes remains elusive. This again under- 
lines that a promising direction is to attempt to find a more 
suitable information-theoretic figure-of-merit that can be used 
in a comparative finite ensemble-size analysis of the different 
QPT schemes. 



3. The Role of Inversion 



It should be noted that in order to reconstruct the unknown 
map £ in a QPT scheme one generally needs to perform an 
inversion operation which here can be understood as Eq. ( fT2b . 
In particular, in the SQPT and AAPT schemes an inversion 
on experimental data is inevitable. This inversion may induce 
an ill-conditioning feature ll25il52ll . i.e., small errors in experi- 
mental outcomes may give rise to large errors in the estimation 
of £, and can sometimes result in non-positive maps. It should 
be stressed that quantum dynamics obtained via the usual pre- 
scription of unitary evolution followed by a partial trace over 
the bath, is always positive when the initial state is a valid 
density matrix. When a positive map is applied outside of 
its positivity domain it will result in non-positive density ma- 
trix. Complete positivity results when in addition one assumes 
a factorized initial system-bath state. Non-complete positiv- 
ity is thus a legitimate feature of correlated initial conditions, 
and non-positivity is a legitimate feature of applying a posi- 
tive map to states outside of its positivity domain |86, 87, 88]. 
The problem with ill-conditioning due to inversion is a of a 
different nature: it is a numerical error that leads to a non- 
positive or non-completely-positive map. This problem, to a 
large extent, can be addressed by supplementary data analysis 
methods, such as ML estimation Ii25i, iSillsl, iMlSSi Is^ dTI] . 
Bayesian state estimation lf58l[59ll60li. an d other reliable regu- 
larization or reconstruction methods I6lll62ll63ll . In principle, 
all known QPT schemes (including DCQD) can be optimized 
by utilizing such statistical error reduction techniques. Here 
we will not delve into the details of such methods, as they are 
applicable on a similar footing to all QPT schemes, and more- 
over, this issue is beyond the scope of the present paper How- 
ever, we would like to emphasize that DCQD is inherently 
more immune against such inversion-amplified errors. The 
diagonal elements of a map, as discussed above, are related in 
DCQD directly to measurement results. For off-diagonal ele- 
ments a large extent of directness also exists. This can easily 
be seen, for example, through the determination of xo3 — see 
Eq. ( fT6] l — in which only the quantities xoo and X33 (already 
obtained from a different experimental configuration) need to 
be used. That is, the formal inversion necessary in DCQD re- 
quires only a small amount of data processing. This, in turn 
implies that inversion-induced errors are amplified less than 
in methods requiring a full inversion. 
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C. Partial Characterization of Quantum Dynamics 

An important and promising advantage of DCQD is for use 
in partial characterization of quantum dynamics, where one 
cannot afford or does not need to carry out a full characteri- 
zation of the quantum system under study, or when one has 
some a priori knowledge about the dynamics. Using indirect 
QPT methods in such situations is generally inefficient, be- 
cause one must apply the whole machinery of the scheme (in- 
cluding its inversion) to obtain the desired partial information 
about the system. On the other hand, the DCQD scheme is 
inherently applicable to the task of partial characterization of 
quantum dynamics. In general, one can substantially reduce 
the total number of measurements when estimating the coher- 
ence elements of the superoperator for only specific subsets 
of the operator basis and/or subsystems of interest. For exam- 
ple, a single ensemble measurement is needed if one wishes to 
identify only the coherence elements xm and xv2 of a partic- 
ular qubit. InRefs. lfl7l[l8ll2lll23ll several examples of partial 
characterization have been demonstrated. For example it was 
demonstrated that DCQD enables the simultaneous determi- 
nation of coarse-grained (semiclassical) physical quantities, 
such as the longitudinal relaxation time Ti and the transver- 
sal relaxation (or dephasing) time T2. Alternative methods for 
efficient selective estimation of quantum dynamical maps has 
been recently developed by utilizing random sampling ll22ll . 
The central idea of these methods is symmetrization of a quan- 
tum channel by randomization, and then efficient estimation 
of gate fidelities. The application of such partial/selective pro- 
cess estimation schemes for efficient Hamiltonian identifica- 
tion of open quantum systems is important per se — ^besides its 
practical implications — and will be presented elsewhere ll23ll . 

VII. CONCLUDING REMARKS 

In the absence of a good, reliable figure-of-merit for the 
performance of QPT schemes, one cannot provide a fully fair 



and decisive comparative analysis. In addition, one should 
also consider the complexity of physical resources associated 
with noisy/imperfect quantum state preparation and measure- 
ments. Nevertheless, in this work we have presented a detailed 
resource-based comparison of ideal quantum process tomog- 
raphy schemes with respect to overall number of experimental 
configurations and elementary quantum operations. In gen- 
eral, SQPT is always the best approach for complete estima- 
tion of quantum dynamical systems when controlled two-body 
interactions are not either available or desirable. However, for 
quantum systems with controllable single- and two-body in- 
teractions, we have shown that the DCQD approach is more 
efficient than SQPT, and all versions of AAPT, in terms of the 
total number of elementary quantum operations required. For 
such systems, DCQD appears attractive for near-term applica- 
tions involving complete verification of small quantum infor- 
mation processing units, especially in trapped-ion and liquid- 
state NMR systems. For example, the number of required ex- 
perimental configurations for systems of 3 or 4 physical qubits 
is reduced from - 5 x 10^ and ~ 6.5 x 10^ in SQPT to 64 and 
256 in DCQD, respectively. Such complete characterization 
of quantum dynamics is essential for verification of quantum 
key distribution procedures, teleportation units (in both quan- 
tum communication and distributed quantum computation), 
quantum repeaters, quantum error-correction procedures, and 
more generally, in any situation in quantum physics where a 
few particles interact amongst themselves and with a common 
environment. 
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